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1. Summary 

Random Matrix Theory has been a unifying approach in physics and math- 
ematics. In these lectures we discuss applications of Random Matrix Theory 
to QCD and emphasize underlying integrable structures. In the first lecture we 
give an overview of QCD, its low-energy limit and the microscopic limit of the 
Dirac spectrum which, as we will see in the second lecture, can be described 
by chiral Random Matrix Theory. The main topic of the third lecture is the 
recent developments on the relation between the QCD partition function and 
integrable hierarchies (in our case the Toda lattice hierarchy). This is an effi- 
cient way to obtain the QCD Dirac spectrum from the low energy limit of the 
QCD partition function. Finally, we will discuss the QCD Dirac spectrum at 
nonzero chemical potential. We will show that the microscopic spectral den- 
sity is given by the replica limit of the Toda lattice equation. Recent results by 
Osborn on the Dirac spectrum of full QCD will be discussed. 
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2. Introduction 

Applications of Random Matrix Theories to problems in physics have a long 
history starting with the idea of Wigner [1] to describe the spacing distribution 
of nuclear levels by an ensemble of real symmetric matrices. Although this is 
the first application of Random Matrix Theory (RMT) to strong interactions, 
applications of RMT to QCD started much later. The first paper that put QCD 
into the context of RMT was the work of 't Hooft on two-dimensional QCD 
in the limit of a large number of colors [2]. It was shown [3] that the combi- 
natorial factors that enter in this large Nc expansion could be obtained from 
matrix integrals. Even today, as we have seen in the lectures by Di Francesco 
[4], this work attracts a great deal of attention. It greatly stimulated the analy- 
sis of a complicated nonlinear theory such as QCD by means of much simpler 
matrix models. Because of the success of the application of RMT to the pla- 
nar expansion of QCD, the hope was that nontrivial results could be derived 
this way. I will mention three well-known results that have emerged from this 
line of thought: the Brezin-Gross-Witten model [5, 6], the Eguchi-Kawai [7] 
reduction, and induced QCD according to Kazakov and Migdal [8]. 

The Wilson loop in lattice QCD without quarks in I-i-l dimensions can be 
reduced to the calculation of the unitary matrix integral 

z{g\N,)=t dUe^-'^(^^^'\ (2.1) 

Juesu{Nc) 

which is known as the Brezin-Gross-Witten model [9, 5, 6]. This reduction 
was generalized to an arbitrary Wilson loop amplitude in the large Nc limit of 
lattice QCD in four dimensions and is known as the Eguchi-Kawai reduction 
[7]. It was shown that Wilson loop amplitudes do not depend on space-time 
and can be obtained from a single plaquette integral. However, this reduction 
is not valid in the weak coupling limit [10]. The idea of induced QCD [8] is 
to induce the plaquette action by a unitary matrix integral. With a vanishing 
Wilson line [11] this approach turned out to be unsuccessful as well. 

The matrix model (2.1) also appears in the low-energy limit of QCD with 
quarks. However, in this case the integral over SU{Nc) is not over the color 
degrees of freedom but rather over the flavor degrees of freedom, and y"^ is 
replaced by mVE/l with S the chiral condensate, m the quark mass and V 
the volume of space time. It coincides with the full QCD partition function in 
a domain where the pion Compton wavelength is much larger than the size of 
the box [12, 13]. In this hmit we have 

^5^.(mSy/2,Ar^) = (:^E^>^n(^^^ + "^)''^)' (2-2) 

where the are the eigenvalues of the Dirac operator. By expanding in powers 
of the inverse mass, one obtains sum rules for the inverse Dirac eigenvalues 
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[13] which put constraints on the Dirac spectrum, but do not determine the 
average spectral density on the scale of the average level spacing (which is 
known as the microscopic spectral density) and other spectral correlators. 

A Random Matrix Theory that describes the fluctuations of the small eigen- 
values of the Dirac operator was introduced in [14, 15]. It was shown that 
chiral RMT is equivalent to the flavor unitary matrix integral (2.1). The spec- 
tral correlation functions were found to be in good agreement with lattice QCD 
simulations (see [16, 17] for a review of lattice results). One argument to un- 
derstand this is that level correlations of complex systems on the scale of the 
average level spacing are universal, i.e. they do not depend on the details of the 
system. This could be shown rigorously in the context of RMT [18]. However, 
it was understood later that the generating function for the Dirac spectrum is 
completely determined by chiral symmetry [19, 20], and its microscopic limit 
does not change if the Dirac operator is replaced by a random matrix. This is 
one of the main topics of these lectures. 

The question has been raised if the quenched spectral density can be ob- 
tained from the limit Nf of (2.2). This procedure is known as the replica 
trick. It has been argued that this limit generally gives the wrong result [21]. 
However, the family of partition functions for different values of A^^ are related 
by the Toda lattice equation [22]. If we take the replica limit of the Toda lattice 
equation [23], or of the corresponding Painlevc equation [24], we obtain the 
correct nonperturbative result. This is a second main topic of these lectures. 

A third main topic of these lectures is the discussion of QCD at nonzero 
baryon chemical potential. In that case the Dirac operator is non-Hermitian 
and its eigenvalues are scattered in the complex plane. We will show that also 
in this case the spectral density can be obtained from the Toda lattice equation 
[25, 26]. Recent analytical results by Osborn [27] for the Dirac spectrum of 
QCD with dynamical quarks at nonzero chemical potential will be discussed. 

We start these lectures with an elementary introduction to QCD and its sym- 
metries. The Dirac spectrum is discussed in section 4. The low energy Umit of 
QCD or partially quenched QCD (see section 5) is equivalent to a RMT with 
the symmetries of QCD that will be introduced in section 6. In section 6 we 
also calculate the microscopic spectral density by means of orthogonal polyno- 
mials and the supersymmetric method. In section 7 we show that this spectral 
density can be obtained from the replica limit of the Toda lattice equation. In 
the same section we connect these results with ideas from the theory of exactly 
solvable systems such as Virasoro constraints, Painleve equations, Backlund 
transformations and the Toda lattice. QCD at nonzero chemical potential is 
discussed in sections 8 and 9. In section 8 we show that the microscopic spec- 
tral density can be obtained from the replica limit of a Toda lattice equation. 
Recent results for full QCD at nonzero chemical potential are discussed in sec- 
tion 9 and concluding remarks are made in section 10. 
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Table 1. Quark Masses 

rriu = 4 McV nic = 1.3 GeV 

md = 8 MeV mb = 4.4 GeV 

= 160 MeV mt = 175 GeV 

Finally, a note about books and reviews on the subject of these lectures. The 
classic RMT text is the book by Mehta [28] which emphasizes the orthogo- 
nal polynomial method. The third edition of this book appeared recently. In 
the book by Forrester [29] the emphasis is on the relation between RMT and 
solvable models and mathematical physics. A comprehensive review of RMT 
is given in [30]. A collection of recent papers on RMT can be found in [31] 
which also contains several reviews. Applications to mesoscopic physics are 
discussed in [32, 33], applications to QCD in [16] and applications to Quan- 
tum Gravity in [34]. Among the pedagogical reviews we mention a general re- 
view of RMT [35] and an introduction to the supersymmetric method [36, 37]. 
These lecture notes overlap in part with lecture notes I prepared for the Latin 
American Summer School in Mexico City [37], where the main emphasis was 
on the supersymmetric method rather than on applications of RMT to QCD. 

3. QCD 

3.1 Introduction 

QCD (Quantum Chromo Dynamics) is the theory of strong interactions that 
describes the world of nucleons, pions and the nuclear force. No experimental 
deviations from this theory have ever been detected. QCD is a theory of quarks 
which interact via gauge bosons known as gluons. In nature we have 6 different 
flavors of quarks which each occur in three colors. Each quark is represented 
by a 4-component Dirac spinor 

ql^, f = l,...,Nf = 6, i = l,...,N, = 3 (3.1) 

with Dirac index /i. In total we have 1 8 quarks (plus an equal number of anti- 
quarks). The gluon fields are represented by the gauge potentials 

Aii, i,j = l,...,iVe = 3, (3.2) 

which, as is the case in electrodynamics, are spin 1 vector fields. The gauge 
fields are Hermitian and traceless; they span the algebra of SU (Nc). In total we 
have 8 massless gluons. The 6 quark flavors are known as up, down, strange, 
charm, bottom and and top. The quark masses are given in Table 1 . Only the 
two lightest quarks are important for low-energy nuclear physics. 

First principle calculations of QCD can be divided into three different groups, 
perturbative QCD, lattice QCD and chiral perturbation theory. The main do- 
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main of applicability of perturbative QCD is for momenta above several GeV. 
Chiral perturbation theory is an expansion in powers of the momentum and 
the pion mass and is only reliable below several hundred MeV. Although lat- 
tice QCD is an exact reformulation of QCD, in practice both the domain of 
low momenta and high momenta cannot be accessed, and its main domain of 
applicability lies somewhere in between the two perturbative schemes. 

The reason that perturbative QCD is applicable at high energies is asymp- 
totic freedom: the coupling constant ^ for momenta p ^ oo. This 
property was instrumental in gaining broad acceptance for QCD as the theory 
of strong interactions and its discovers were awarded this years Nobel prize. 

A second important property of QCD is confinement, meaning that only 
color singlets can exist as asymptotic states. This empirically known property 
has been confirmed by lattice QCD calculations. However, a first principle 
proof of the existence of a mass gap in QCD is still lacking even in the absence 
of quarks. Because of confinement the lightest particles of the theory are not 
quarks or gluons but rather composite mesons. 

This brings us to the third important property of QCD: chiral symmetry. At 
low temperatures and density chiral symmetry is broken spontaneously which, 
according to Goldstone's, theorem gives rises to massless Goldstone bosons. 
Because the chiral symmetry is slightly broken by the light quark masses, the 
Goldstone bosons are not exactly massless, but the mass of the pions of 135- 
138 MeV is an order of magnitude less than a typical QCD scale of about 1 
GeV. This justifies a systematic expansion in the pion mass and the momenta 
known as chiral perturbation theory. 

A fourth important property of QCD is that a first principle nonperturba- 
tive lattice formulation can be simulated numerically. This allows us to com- 
pute nonperturbative observables such as for example the nucleon mass and 
/?-meson mass. Without lattice QCD we would have had only a small number 
of first principle nonperturbative results and the validity of QCD in this domain 
would still have been a big question mark. 

3.2 The QCD partition function 

The QCD partition function in a box of volume V3 = can be expressed 
in terms of the eigenvalues of the QCD Hamiltonian as 



where (3 is the inverse temperature. At low temperatures, (/3 — > cxd), the parti- 
tion function is dominated by the lightest states of the theory, namely the vac- 
uum state, with an energy density of Eq/Vs and massless excitations thereof. 
The partition function ^Q^D rewritten as a Euclidean functional inte- 




(3.3) 



k 
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gral over the nonabelian gauge fields A^, 

N 



JV f 

Z^^^{M) = jdA^ll det(L> + m/)e-^™, (3.4) 



/=i 

where is the Yang-Mills action given by 



S^"" = I dH^F^:.' - (3-5) 



The field strength and its dual are given by 



(3.6) 



The fabc the structure constants of the gauge group SU{Nc). The gauge 
fields are denoted by = A^a'^, where are the generators of the gauge 
group. The integral v = J d^xF^^^FI^^ is a topological invariant, i.e. it 
does not change under continuous transformations of the gauge fields. An 
important class of field configurations are instantons. These are topological 
nontrivial field configurations that minimize the classical action. They are clas- 
sified according to their topological charge u. The parameter 9 is known as the 
0-angle. Experimentally, its value is consistent with zero. In (3.4), the mass 
matrix is diagonal, M = diag(mi, • • • , niNf), but below we will also consider 
a general mass matrix. The anti-Hermitian Dirac operator in (3.4) is given by 

D = j^{d^ + iA^), (3.7) 

where the are the Euchdean Dirac matrices with anti-commutation relation 
{in, 7i/} = 2(5^i^. In the chiral representation the 7-matrices are given by 

In this representation the Dirac operator has the structure 

id 
it 



The integration measure is defined by discretizing space-time 

dA''^ = l[dA''^{x). (3.10) 



A particular popular discretization is the lattice discretization where the QCD 
action is discretized on a hyper-cubic lattice with spacing a. The discussion of 
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lattice QCD would be a lecture by itself. For the interested reader we recom- 
mended several excellent textbooks on the subject [38, 39, 40]. 

A field theory is obtained by taking the continuum limit, i.e. the limit of zero 
lattice spacing a for the integration measure discussed above. This limit only 
exists if we simultaneously adjust the coupling constant, i.e. g — ^ g{a). If such 
limit exists the field theory is called renormalizable. For QCD g{a) approaches 
zero in the continuum limit, a property known as asymptotic freedom. 

We will be mainly interested in the eigenvalues of the Dirac operator and 
how they fluctuate for gauge fields distributed according to the QCD action. 
We will show that below a well-defined scale the fluctuations of the Dirac 
eigenvalues are given by a RMT with the global symmetries of the QCD. 

3.3 Symmetries of the QCD Partition Function 

It is well-known that the QCD action is greatly constrained by gauge sym- 
metry, Poincarc invariance and renormalizability. These symmetries determine 
the structure of the Dirac operator and are essential for its infrared spectral 
properties. In this section we will discuss the global symmetries of the Eu- 
clidean Dirac operator In particular, the chiral symmetry, the flavor symmetry 
and the anti-unitary symmetry of the continuum Dirac operator are discussed. 

3.3.1 Axial Symmetry. The axial symmetry, or the C/a(1) symmetry, 
can be expressed as the anti-commutation relation 



This implies that all nonzero eigenvalues occur in pairs itiAfc with eigenfunc- 
tions given by (pk and ^^4>k- If = the possibility exists that 75(/)fe ~ so 
that Afc = is an unpaired eigenvalue. According to the Atiyah-Singer index 
theorem, the total number of such zero eigenvalues is a topological invariant, 
i.e., it does not change under continuous transformations of the gauge field 
configuration. Indeed, this possibility is realized by the field of an instanton 
which is a solution of the classical equations of motion. On the other hand, it 
cannot be excluded that A^ = while (t>k and 75 (/>a; are hnearly independent. 
However, this imposes additional constraints on the gauge fields that will be 
violated by infinitesimal deformations. Generically, such situation does not 
occur. 

In a decomposition according to the total number of topological zero modes, 
the QCD partition function can be written as 



{75,1?} = 0. 



(3.11) 



00 




(3.12) 
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where 

ZQCD(M) = ([[m}l[i>^l + m})),. (3.13) 

/ k 

Here, {■ ■ ■)u denotes the average over gauge-field configurations with topo- 
logical charge u weighted by the Yang-Mills action. If we introduce right- 
handed and left-handed masses as complex conjugated masses we find that the 
9 dependence of the QCD partition function is only through the combination 
777,e'^/^/. This property can be used to obtain the 6'-dependence of the low- 
energy effective partition function. 

3.3.2 Flavor Symmetry. A second important global symmetry is the 
flavor symmetry. This symmetry can be best explained by writing the fermion 
determinant in the QCD partition function as a functional integral over Grass- 
mann variables, 

ndet(L> + mf)= [ dV'dVie/'^'^S/=i^'(-°+"^^)^'. (3.14) 

In a chiral basis with i/jr = ^^i/^r and f/'L = —Ibi^L, the exponent can be 
rewritten as 

/=i 

(3.15) 

To better illuminate the transformation properties of the partition function 
we have replaced the diagonal mass matrix by Mrl and Mlr. 
For m/ = we have the symmetry 

IpL^ULlpL, i'L^i'lUl^, ''JjR^URtpR, tpR ^ iPrU]^^. (3.16) 

The only condition to be imposed on U and V is that their inverse exists. If the 
number of left-handed modes is equal to the number of right-handed modes 
we thus have an invariance under GlR{Nf) x GlL{Nf), where Gl{Nf) is the 
group of complex Nj x Nj matrices with nonzero determinant. However, if 
the number of left-handed modes is not equal to the number of right-handed 
modes, the axial-symmetry group is broken to an Sl{Nf) subgroup whereas 
the vector symmetry with Ul = Ur remains unbroken. For m/ = the flavor 
symmetry is thus broken explicitly to Glv{Nf) x SlA{Nf) by instantons or the 
anomaly. A subgroup of Glv{Nf) corresponds tobaryon number con- 

servation and is usually not considered when flavor symmetries are discussed. 
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What is much more important, though, is the spontaneous breaking of the 
axial flavor symmetry. From lattice QCD simulations and phenomenological 
arguments we know that the expectation value (V'V') = (V'-rV'l) + (^lV'r) ~ 
— (240 MeV)^ in the vacuum state of QCD instead of the symmetric possibihty 
{'(pip) = 0. Phenomenologically, this is known because the pions are much 
lighter than the a mesons. The spontaneous breaking of the axial symmetry 
also follows from the absence of parity doublets. For example, the pion mass 
and the oq mass are very different (mTr = 135 MeV and rrij = 980 MeV). 

For fermionic quarks there is no need to extend the symmetry group to 
Glji{Nf) X GlL{Nf). In that case, the flavor symmetry is spontaneously bro- 
ken according to SUR{Nf) x SUiiNf) SUv{Nf). For bosonic quarks, 
which enter in the generating function for the Dirac spectrum, it will be shown 
in the next section that it is essential to consider the complex extension of 
SU(Nf ). Notice that the complex extension of a symmetry group does not 
change the number of conserved currents. 

On easily verifies that {'ij^'tp) is only invariant for Ul = Ur. The vacuum 
state thus breaks the chiral symmetry down to Gl\/{Nf). In agreement with 
the Vafa-Witten theorem [41] only the axial symmetries can be broken spon- 
taneously. We also observe that the complete axial group is broken which is 
known as the maximum breaking [42] of chiral symmetry. 

3.3.3 Flavor Symmetry for Bosonic Quarks. For bosonic quarks the 
Goldstone bosons cannot be parameterized by a unitary matrix. The reason is 

that symmetry transformations have to be consistent with the convergence of 
the bosonic integrals. Let us consider the case of one bosonic flavor. Then 

(3.17) 

so that the integral is convergent for Re(m) > 0. The most general flavor 
symmetry group of the action in (3.17) is Gl{2) that can be parameterized as 

U = e"V with = H and VV'^ = 1. (3.18) 

For [/ to be a symmetry transformation for m = we require that 




(3.19) 



so that H has to be a multiple of 0-3, and V has to be a multiple of the identity. 
The transformations V in (3.18) are not broken by the mass term and therefore 
represent the vector symmetry. Only the symmetry transformation exp(scr3) is 
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broken by the mass term so that the axial transformations are parameterized by 
U = (^^^ J^s^ with sG (-00,00). (3.20) 
For Nf bosonic flavors the axial transformations are parameterized by 

U= (^^^ ^ j with = H, (3.21) 

which is the coset Gl{Nf)/U{Nf). 

3.4 Anti-Unitary Symmetries and the Three-fold Way 

The QCD partition function with three or more colors in the fundamental 
representations has no anti-unitary symmetries. As will be discussed below, 
for two colors with fundamental fermions and for adjoint fermions, the Dirac 
operator has an anti-unitary symmetry. The classification of the QCD Dirac 
operator according to anti-unitary symmetries was introduced in [15]. 

3.4.1 Anti-Unitary symmetries and the Dyson index. The value of 

the Dyson index is determined by the anti-unitary symmetries of the system. 
If there are no anti-unitary symmetries the Hamiltonian is Hermitian and the 
value of the Dyson index is Pd = 2. 

An anti-unitary symmetry operator, which can always be written as A = 
U K with U unitary and K the complex conjugation operator, commutes with 
the Hamiltonian of the system 

[H, UK] = 0. (3.22) 

We can distinguish two possibilities, 

{UKf = 1 or {UKf = -1, (3.23) 

corresponding to /3d = 1 and [3d = 4, respectively. The argument goes as 
follows. The symmetry operator = {UK)'^ = UU* is unitary, and in an 
irreducible subspace, it is necessarily a multiple of the identity, UU* = Al. 
Because of this relation, U and U* commute so that A is real. By unitarity we 
have |A| = 1 which yields A = ±1. 

When = 1 it is always possible to find a basis in which the Hamiltonian 
is real. Starting with basis vector 0i we can construct ipi = (j)i + U K(f)i. Then 
choose (f)2 perpendicular to V'l and define '4>2 = (f>2 + UK(f)2 with 

{(t>2 + UKct>2, ^1) = {UKct>2, ^1) = {{UKfct>2, UK^iY = {ct>2, V'l)* = 0. 
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The next basis vector is found by choosing ^3 perpendicular to tpi and 
etc. . In this basis the Hamiltonian is real 

Hki = i^k, HiJi) = (UKi^k, UKH^i)* = {iPk, HUK^iY = {^k, H^lT 
= mi- (3.24) 

The best known anti-unitary operator in this class is the time-reversal oper- 
ator for which U is the identity matrix. 

In the case {UK)"^ = —1 all eigenvalues of the Hamiltonian are doubly 
degenerate. This can be shown as follows. If 0^ is and eigenvector with eigen- 
value Afc, then it follows from (3.22) that also UKcf)^ is an eigenvector of the 
Hamiltonian with the same eigenvalue. The important thing is that this eigen- 
vector is perpendicular to 0fc [13], 

((^fe, UK<t>k) = {UK(l>k, {UKfcPkr = -(<^fe, UK^k)- (3.25) 

In this case it is possible to construct a basis for which the Hamiltonian ma- 
trix can be organized into real quaternions [43]. The eigenvalues of a Hermi- 
tian quaternion real matrix are quaternion scalars, and the eigenvalues of the 
original matrix are thus doubly degenerate in agreement with (3.25). The best 
known example in this class is the Kramers degeneracy for time reversal invari- 
ant systems with half-integer spin but no rotational invariance. For example, 
for spin | the time reversal operator is given by with {Ka2)^ = —1. 
Next we will discuss the anti-unitary symmetries of the QCD Dirac operator. 

3.4.2 QCD in the Fundamental Representation. For three or more 
colors, QCD in the fundamental representation does not have any anti-unitary 
symmetries and fdo = 2. QCD with two colors is exceptional. The reason is 
the pseudo-reality of SU{2): 

A; = {j:A'^^^r = -r2A,r2, (3.26) 
k ^ 

where the r/j are the Pauli matrices acting in color space. From the explicit 
representation for the 7-matrices it follows that 

7/^ = 72747/^7274- (3.27) 

For the Dirac operator iD = ij^dfj, + j^Af^ we thus have 

[KC^5r2,D]=0, (3.28) 

where K is the complex conjugation operator and C = 7274 is the charge 
conjugation matrix. Because (^^6*75X2)^ = 1 we have that ftjy = 1. Using the 
argument of Eq. (3.24) a basis can be constructed such that the Dirac matrix is 
real for any A^. 
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3.4.3 QCD in the Adjoint Representation. For QCD with gauge fields 
in the adjoint representation the Dirac operator is given by 

D = -ff^d^Shc + r'^'l^Aa^, (3.29) 

where the /"'"^ denote the structure constants of the gauge group. Because of 
the complex conjugation property of the 7-matrices we have that 

b27a5K,D]=0. (3.30) 

One easily verifies that in this case 

(727475^)' = -1, (3.31) 

so that the eigenvalues of D axe doubly degenerate (see section 3.4.1). This 
corresponds to the case /Si? = 4, so that it is possible to organize the matrix 
elements of the Dirac operator into real quaternions. 

4. The Dirac Spectrum in QCD 

In this section we show that the smallest eigenvalues of the QCD Dirac 
operator are related to the chiral condensate by means of the Banks-Casher 
relation. This result is used to define the microscopic spectral density. 

4.1 Banks-Casher Relation 

The order parameter of the chiral phase transition, {'ijjip), is nonzero only 
below a critical temperature. As was shown by Banks and Casher [44], {'ijjtp) 
is directly related to the eigenvalue density of the QCD Dirac operator per unit 
four-volume 

S=|(V;V)|=lim^^^. (4.1) 
For eigenvalues {Afc} the average spectral density is given by 

p{X) = {Y,Si\-Xk)). (4.2) 

k 

To show the Banks-Casher relation we study the resolvent defined by 

It can be interpreted as the electric field at z of charges at iXk- Using this 
analogy it is clear that the resolvent changes sign if z crosses the imaginary 
axis. Let us look at this in more detail. A typical Dirac spectrum is shown in 
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Figure 1. 

A typical Dirac spectrum. To derive the Banks-Casher relation we integrate the resolvent over 
the rectangular contour in this figure. (Figure taken from [45].) 



Fig. 1 . The average number of eigenvalues in the rectangular contour in this 
figure is p{\)l. If we integrate the resolvent along this contour we find 

j> G{z) = il{G{i\ + e) - G{i\ - e)) = 27rip(A)/, (4.4) 

where the second identity follows from Cauchy's theorem. Using the symme- 
try of the spectrum we obtain for e — > 

ReG(iA + e) = 7rp(A). (4.5) 

Near the center of the spectrum the imaginary part of the resolvent is negUgi- 
ble. Using that the chiral condensate is related to the resolvent by 

(^/S-^) = - lim lim ^G{m), (4.6) 

immediately results in the Banks-Casher relation (4.1). The order of the limits 
in (4.1) is important. First we take the thermodynamic Umit, next the chiral 
limit (and, finally, the continuum limit). 

The resolvent of the QCD Dirac spectrum can be obtained from 

G{z,z',mf)= ^ logZPi(z,z',m/), (4.7) 

az 2=^/ 

with the so called partially quenched QCD partition function given by 

Zl\z,z\mj)= j dA ^ + n det(i^ + m,) e'^"" . (4.8) 
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For z = z' this partition function coincides with the QCD partition function. 

In addition to the regular quarks, the partition function (4.8) has additional 
bosonic and fermionic ghost quarks. Our aim is to find the chiral Lagrangian 
corresponding to (4.8). If we are successful, we have succeeded in deriving a 
generating function for the infrared Umit of the QCD Dirac spectrum. 

4.2 Microscopic Spectral Density 

An important consequence of the Bank-Casher formula (4.1) is that the 
eigenvalues near zero virtuality are spaced as 

AA = l/p(0) = vr/El^. (4.9) 

For the average position of the smallest nonzero eigenvalue we obtain the esti- 
mate 

Amin = Tt/S^. (4.10) 

This should be contrasted with the eigenvalue spectrum of the non-interacting 
Dirac operator. Then the eigenvalues are those of a free Dirac particle in a 
box with eigenvalue spacing equal to AA ~ for the eigenvalues near 

A = 0. Clearly, the presence of gauge fields leads to a strong modification of 
the spectrum near zero virtuality. Strong interactions result in the coupling of 
many degrees of freedom leading to extended states and correlated eigenvalues. 
Because of asymptotic freedom, the spectral density of the Dirac operator for 
large A behaves as V)? . In Fig. 2 we show a plot of a typical average spectral 
density of the QCD Dirac operator for A > 0. The spectral density for negative 
A is obtained by reflection with respect to the y-axis. More discussion of this 
figure will be given in section 5.3. 

Because the eigenvalues near zero are spaced as ~ it is natural to 

introduce the microscopic spectral density [14] 

1 u 

p,(u) = lim — p(— ) with u = \VT.. (4.11) 

We expect that this limit exists and converges to a universal function which is 
determined by the global symmetries of the QCD Dirac operator. In section 6, 
we will calculate ps (u) both for the simplest theory in this universality class, 
which is chiral Random Matrix Theory (chRMT), and for the partial quenched 
chiral Lagrangian which describes the low-energy Umit of the QCD partition 
function. We will find that the two results coincide below the Thouless energy. 

5. Low Energy Limit of QCD 

In this section we derive the chiral Lagrangian that provides an exact de- 
scription of QCD at low energies. 
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Figure 2. 

Schematic picture of the average spectral density of QCD Dirac operator. (Taken from [46].) 

5.1 The chiral Lagrangian 

For light quarks the low energy limit of QCD is well understood. Because of 
confinement, the only light degrees of freedom are the pseudo-scalar mesons, 
that interact according to a chiral Lagrangian. To lowest order in the quark 
masses and the momenta, the chiral Lagrangian is completely determined by 
chiral symmetry and Lorentz invariance. In the case of Nf light quarks with 
chiral symmetry breaking according to SUiiNf) x SUR{Nf) SUv{Nf) 
the Goldstone fields (pseudo-scalar mesons) are given by C/ G SU (Nf). Under 
an SUiiNf) x SUR{Nf) transformation of the quark fields given in (3.16), 
the Goldstone fields U transform in the same way as the chiral condensate 

U^UrUUi\ (5.1) 

The symmetry (3.16) is broken the mass term. However, the full symmetry can 
be restored if we also transform the mass term as 

Mrl ^ UrMrlUZ\ Mlr ^ UlMlrU^K (5.2) 

The low energy effective theory should have the same invariance properties. 
To leading order it contains terms of second order in the momenta and of first 
order in the quark mass matrix. The invariant terms are: 

TY{d^Ud^U^), lh:{MRLU^), Tr{MLRU). (5.3) 

Since the QCD partition function is invariant under Mrl ^ Mlr, the effec- 
tive partition function should also have this symmetry. The action of the Gold- 
stone fields is therefore given by the so called Weinberg Lagrangian [47, 48] 

Ces{U) = ^ Ti^d^Ud^U"^) - I T:v{MrlU^ + MlrU), (5.4) 
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where F is the pion decay constant, and S is the chiral condensate. The Gold- 
stone fields can be parametrized as [/ = cxp(i V2n„t"/F), with the genera- 
tors of SU{Nf) normahzed as Tr t"-t^ = 6"-^. This chiral Lagrangian has been 
used extensively for the analysis of pion-pion scattering amplitudes [48]. 
To lowest order in the pion fields we find for equal quark masses m 

CesiU) = ^d^U'^d'^U'' - NfT.m + ^n»n». (5.5) 

This results in the pion propagator l/(p^ -|- m^) with pion mass given by the 
Gellmann-Oakes-Renner relation 

ml = -y^. (5.6) 

It also illustrates the identification of S as the chiral condensate. 

5.2 The Low Energy Limit of Zp'i 

The low-energy limit of the partially quenched QCD partition function can 

be derived along the same lines as the derivation of the chiral Lagrangian ob- 
tained in previous section. In this case, ignoring convergence questions for the 
moment, the global flavor symmetry of (4.8) is given by by the supergroup 

GlRiNf + l\l)xGlL{Nf + l\l). (5.7) 

This reflects that we have Nf + 1 fermionic quarks and 1 bosonic quark. We 
already have seen that convergence requirements restrict the axial symmetty 
for bosonic quarks to GZ (!)/[/(!). Although the axial flavor symmetry group 
of the fermionic quarks is not a priori determined by convergence requirements 
we will see in this section that supersymmetry necessarily imposes that this 
symmetry group is compact, i.e. equal to U {Nf). 

Under transformation (5.7) the quarks fields with Nf + 1 fermionic compo- 
nents and one bosonic component, transform as (Ur, Ul G Gl{Nf -I- 1|1)) 

ipR Uri/jr, tpL UliPl, ipR iPrUr^, ipL iPlU£^. (5.8) 

The subscripts refer to the right-handed (R) or left-handed (L) quarks. For 
M = and z/ = this is a symmetry of the QCD action. For M / this 
symmetry can be restored if we also transform the mass term according to 

Mrl ^ UrMrlUI\ Mlr ^ UlMlrU^\ (5.9) 

In the sector of topological charge v the partially quenched partition function 
transforms as 



ZI\Mrl, Mlr) ^ ^dei^'pRUL]-^ Z^^Mrl, Mlr). (5.10) 
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The Goldstone bosons corresponding to the breaking of the axial subgroup 
GlA{Nf + 1|1) transform as Q -> UrQU£^. If we factorize the Goldstone 
fields into the zero momentum modes Qo and the nonzero momentum modes 

Q{x) as 

Q = QoQ{x), (5.11) 

one can easily show that the low energy effective partition function with the 
above transformation properties is given by 

ZPi(M)= / dQSdet''{Qo)e-f'^^''^'"'^^\ (5.12) 

JQeGliNf+l\l) 

where (see [32, 37] for definitions of the superdeterminant Sdet and and the 
supertrace Str) 

^""""{Q) = ^Strd^Q-^d^Q - |str(MRiQ-i) - |str(MiRQ). (5.13) 

We already have seen that the boson-boson block of Gl{Nf+l\l) is (1)/C/(1) 
If we parameterize the field Q as 

Q = e^k'^'^'"^/^, (5.14) 

with Tjfc the generators of G{Nf + 1| 1), to second order in the Goldstone fields 
the mass term is given by Str{Y,MJ2k'^k'^k/^'^)- Let us take M diagonal 
positive definite. Because of the supertrace there is a relative minus sign be- 
tween the boson-boson and fermion-fermion modes. The boson-boson modes 
are noncompact and require that the overall minus sign of the mass term is neg- 
ative. In order to avoid tachyonic fermion-fermion Goldstone modes, we have 
to compensate the minus sign of the supertrace. This can be done by choos- 
ing the parameters that multiply the fermion-fermion generators purely imag- 
inary. This corresponds to a compact parametrization of the fermion-fermion 
Goldstone manifold. This integration manifold is the maximum Riemannian 
submanifold [49] of Gl{Nf + 1|1) and will be denoted by Gl{Nf + 1|1). 

5.3 The Mesoscopic Limit of QCD 

In chiral perturbation theory, the different domains of validity where ana- 
lyzed by Gasser and Leutwyler [12]. A similar analysis applies to partially 
quenched chiral perturbation theory [51]. The idea is as follows. The Q field 
can be decomposed as [12] 



(5.15) 
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where Qq is a constant (zero-momentum) field. For momenta p = nk/L with 
k integer, the kinetic term of the ^ fields behaves as 

^di,tjj''{x)di,tjj''{x) ~ L-V(a;)V'"(2:). (5.16) 

We observe that the magnitude of the fluctuations of the ip°'{x) fields are of 
order 1/L which justifies a perturbative expansion of exp(i^(a;)). The fluctu- 
ations of the zero momentum modes are only limited by the mass term 

^FEStrM(Qo + Qo')- (5-17) 

For quark masses m S> 1 /VTf, the field Qo fluctuates close to the identity and 
the Qo field can be expanded around the identity as well. If m <C Aqcd we 
are in the domain of chiral perturbation theory. For 

the fluctuations of the zero momentum modes dominate the fluctuations of the 
nonzero momentum modes, and only the contribution from the zero momen- 
tum modes has to be taken into account for the calculation of an observable. In 
this hmit the so called finite volume partition function is given by [12, 13] 

rM,0)~ / d^7e^^I^^^^^'^^''"^ (5.19) 

^ Ju<=SU{Nf) 

where the ^-dependence follows from the dependence of the QCD partition 
function on the combination me^^^^f only (see section 3.3.1). We emphasize 
that any theory with the same pattern of chiral symmetry breaking as QCD can 
be reduced to the same extreme infrared limit. 

The effective partition function at fixed u follows by Fourier inversion [13] 

1 /■^''' 

ZfNAM) = T^ dee-'''^Z''^{M,e). (5.20) 
' ' 27r Jo 

Combining the integral over SU (Nf) and the integral over U{1) we find that 
Zf^(M)= [ det'^(C/)e^^^^^^^^ (5.21) 

' ' Ju{Nf) 

The same arguments apply to the partially quenched chiral Lagrangian. 
There is an important difference. The mass of the ghost-quarks is an exter- 
nal parameter which can take on any value we wish. Indeed, the mass of the 
Goldstone modes containing these quarks is given by 

Ml = (5.22) 
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Therefore, independent of the quark masses there is always a domain where 
the fluctuations of the zero momentum modes dominate the fluctuations of the 
nonzero momentum modes. This domain is given by [51] 

z < = mc. (5.23) 

In this domain, the Compton wavelength of the Goldstone bosons with mass 
Mzz is much larger than the size of the box. Because the time scale conjugate 
to Txic is of the order of the diffusion time across the length of the box, this 
domain is known as the ergodic domain. For the non-Goldstone modes not to 
contribute to the partition function we have to require that L » 1/Aqcd- 
In the Dirac spectrum we can thus distinguish three important scales: 

Amin < mc < Aqcd (5.24) 

(see Fig. 2). For z <C nic we are in the zero momentum sector of the theory. 
If z is of the order of Amin or less we have to take into account quantum fluc- 
tuations to all orders. For Amin ^ z ^ rric, the integral over zero modes can 
be calculated perturbatively by a loop expansion. For mc <C z <C Aqcd, chi- 
ral perturbation theory still applies, but the zero momentum modes no longer 
dominate the partition function. For z 2> Aqcd, chiral perturbation theory is 
not appUcable to the spectrum of the Dirac operator. 

In the ergodic domain the QCD partition function in the sector of topological 
charge u is given by [19] 

Zl^{M)=[ dQSdet'^Qe^?str(MQ+MQ-i)_ ^^^5) 

JQ&Gl{Ni+l\l) 

The number of QCD Dirac eigenvalues that is described by this partition func- 
tion is of the order mc/AA = F^L'^. This number increases linearly in Nc for 
00 which was recently found in lattice simulations [50]. 
In section (6.3.2) we will study this partition function in the quenched limit 
(A^^ = 0) and show it coincides with the chRMT result [19, 20]. 

5.3.1 Comparison to Disordered Systems. In the book by Efetov [32] 
it is shown that the diffusion of electrons in a disordered medium can be de- 
scribed by the effective action 

F{Q) = J d'^x[^DTr{VQf - ™TVAQ], (5.26) 

where Q are the Goldstone fields, v is the density of states, D is the diffusion 
constant and oj is the energy difference between the advanced and the retarded 
Green's functions. The matrix A is a diagonal matrix with matrix elements ±1 



25 



corresponding to the causal character of the Green's functions. The Goldstone 
bosons arise because of the spontaneous breaking of the symmetry between 
the advanced and retarded Green's functions. 

If we compare tliis effective action to the chiral Lagrangian (5.13) we can 
make the identification 



TTZ/D TTUJV MS p{E) 



4 8 ' 4 2 ' V 

which can be rewritten as 



(5.27) 



M ^ -, S^TTZ/, ^ 2~' ^^''^^^ 

The domain where the kinetic term factorizes from the partition function is 
therefore given by [51, 52] 



« — - ^ -. (5.29) 

In the theory of disordered mesoscopic systems the corresponding energy scale 
is known as the Thouless energy. It is defined by [53, 54]: 

UD 

Ec = jj, (5.30) 

The time conjugate to Ec is the time scale over which an electron diffuses 
across the sample. Therefore, the domain where hu ^ Ec is known as the 
ergodic domain. The time scale in mesoscopic physics corresponding to IVqcd 
is the elastic scattering time Tg. The domain in between Ec and h/Teis known 
as the diffusive domain. This domain is characterized by diffusive motion of 
electrons in the disordered sample described by the Lagrangian (5.26). 

6. Chiral RMT and the QCD Dirac Spectrum 

6.1 The chiral ensembles 

The chiral ensembles are defined as the ensembles of N x N Hermitian 
matrices with block structure [14, 15] 

and probabiUty distribution given by (for equal quark masses m) 

P{C)dC = Afdet^f {D + m)e-^'^^^^dC. (6.2) 

The integration measure dC is the product of differentials of the independent 
parts of the matrix elements of C, and Nf is a real parameter (corresponding to 



F2 D 
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the number of quark flavors in QCD). The matrix C is a rectangular nx{n+iy) 
matrix and N = 2n + u. The nonzero eigenvalues of the matrix D occur in 
pairs ±Afe. This can be seen as follows. If 

n{l)^^{:) .he„ -J.-a(_«J. ,6.3) 

Generically, the matrix D in (6.1) has exactly \u\ zero eigenvalues. For this 

reason, u is identified as the topological quantum number. The normalization 
constant of the probabihty distribution is denoted by Af. We can distinguish 
ensembles with real, complex, or quaternion real matrix elements. They are 
denoted by /3d = ^, Pd = 2, and pjj = 4, respectively. In addition to 
the global symmetries of QCD, this partition function has the large unitary 
invariance 

C^UCV-^, (6.4) 

where U and V are orthogonal, unitary, or symplectic matrices, respectively. 
Therefore, the corresponding ensembles are known as the chiral Gaussian Or- 
thogonal Ensemble (chGOE), the chiral Gaussian Unitary Ensemble (chGUE), 
and the chiral Gaussian Symplectic Ensemble (chGSE), in this order. 
Using the invariance (6.4) it is always possible to decompose C as 

C = C/AF"S (6.5) 

where A is a diagonal matrix with > 0. The joint probability distribution for 
the eigenvalues is obtained by transforming to A, U and V as new integration 
variables. The Jacobian is given by 

k k<l 
resulting in the joint eigenvalue distribution 

P{{X})d{\} = Af\A{{X^})f'^ H X^{Xl + m2)^/e-^^°^^/^dAfc, (6.7) 

k 

where a = /3u — 1 + /^iji^. We note that the distribution of the eigenvectors 
factorizes from the distribution of the eigenvalues. 

6.2 Mathematical Methods 

In this section we will discuss the orthogonal polynomial method, the re- 
solvent expansion method, the replica trick and the supersymmetric method. 
These methods are widely used in Random Matrix Theory. 
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6.2.1 Resolvent Expansion Methods. These methods are based on ex- 
panding the resolvent in a geometric series 

G{z) = (TV^-) =N- + {Tr-H-) + {Tr-H-H-) + ■■■ . (6.8) 
z — H z z z z z z 

In the large N hmit the averages are given by a sum of planar diagrams. Let 
us illustrate this for the GUE. In this case the "propagator" is given by 



{HijHki) — —Siidjk- 



(6.9) 



For example, as was explained in the course of Di Francesco [4], for TrH^ 
term we have two planar diagrams of order and one diagram of order N'^. 

6.2.2 The Orthogonal Polynomial Method. The oldest method is the 
orthogonal polynomial method [28]. In principle, one obtains expressions that 

are exact for finite size matrices. The drawback of this method is that it requires 
a probability distribution that is invariant under basis change of the random 
matrix. In general the probabiUty density can be written as 



P{X1, ■■■,Xn) = A/^°({xfe}) n w{Xk), 

k=l 



(6.10) 



where w{x) is a weight function and the Vandermonde determinant is given by 

A({xfc}) = n(^fc-^0- (6.11) 

k>l 



The method is based on the identity 



A(K}) 



1 

Xl 



Xl 



x: 



1 

n-1 



Poixi) 
Pi{xi) 



n-1 



(xi) 



PoiXn) 
Pl{Xn) 

Pn—liXn) 



where the Pyt are monic orthogonal polynomials defined by 
dxw{x)Pk{x)Pi{x) = hkSki- 



,(6.12) 



(6.13) 



Because of these relations, integrals over the eigenvalues can be performed by 
means of orthogonality relations. In the next section we illustrate this method 
by the calculation of the microscopic spectral density for the chGUE. 
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6.2.3 The Replica Trick. The replica trick is based on the identity 

G{z) = i(TV^) = limi^5,(det'^(ii; + .)). (6.14) 

The recipe is to calculate the partition function for positive or negative integer 
values of r and then analytically continue to r = 0. For positive (negative) 
integer values of r the average determinant can be calculated by rewriting it 
as a Grassmann (complex) Gaussian integral. Then D appears Unear in the 
exponent which allows us to perform the average for a Gaussian distribution of 
D. The replica trick works without problems for perturbative calculations [55] 
but usually fails in nonperturbative calculations [21]. As example consider the 
following expression for the the modified Bessel function luiz): 

1 sin zvTT /"^ 
I^{z) = - e'^^^'^cosuOde-- -/ e-^"°^^*-^*dt. (6.15) 

TT io Jo 

We would have missed the second term if we calculate the Bessel function only 
for integer values of u. The replica trick can be made to work if we consider 
a family of partition functions related by a Toda lattice equation. This will be 
discussed in detail in the next two lectures. 

6.2.4 The SupersymmetriclMethod. The supersymmetric method [76] 
is based on the identity 

(6.16) 

z'=z 

The determinant can be written as a Grassmann integral and the inverse deter- 
minant as a complex integral. For z' = z this partition function has an exact 
supersymmetry. The advantage of this method is that it is mathematically rig- 
orous, but it requires a deep understanding of super mathematics. For example, 
finite expressions can be obtained from singular terms that do not depend on 
the Grassmann variables (and are zero upon integration). For a discussion of 
these so-called Efetov-Wegner terms we refer to the hterature [56, 37]. 

6.3 The IMicroscopic Spectral Density of the chGUE 

In this subsection we calculate the microscopic spectral density by the or- 
thogonal polynomial method [57] and the supersymmetric method [19, 20]. 

6.3.1 Orthogonal Polynomials. For the chGUE the joint probability 
distribution only depends on the square of the eigenvalues. In terms of Xk = 
the weight function for m/ = is given by 



1 /det^W + z)_\ 



w{x) = (a;nS^)"e ' 



(6.17) 
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with a = Nf + The monic orthogonal polynomials corresponding to this 
weight function can be expressed in terms Laguerre polynomials 



(6.18) 



with normalization constants hk given by = k\{k + a)!/(nS^)^'^+^. The 
eigenvalue density is given by (with c a constant) 



/n 
JJ [w{xk)dxk 

k=2 



Po{Xn) 



P„_l(xi) ••• Pn-liXn) 

n 

n— 1) (-^n) 



k=2 



X^7r(0) (•^l)---^7r(n-l)(a^n), 



n— 1 n— 1 



{n-\)\W}HY.-Pl{x^)w{x{) 



(6.19) 



The microscopic spectral density is obtained by taking the limit n ^ cxd for 
fixed z = 2nSA = 2nT,^/x. In this limit the weight function is given by 
w{x) = (xS^n)", and the Laguerre polynomials behave as 



Ll{xT,'^n) A;"(xS2rafe)-"/2ja(2S\/^). 



(6.20) 



In the limit n ^ oo, the sum can be replaced by an integral resulting in the 
microscopic spectral density [57] 



p{z) ~ zJ^Ja 

k=0 

= 2Z{JI{Z) - Ja+l{z)Ja-l{z)) 




tdtJlizt), 



(6.21) 



6.3.2 Supersymmetric Method. In this section we evaluate the re- 
solvent of QCD for the simplest case of Nf = and = in the domain 
z ^ F'^ jYuI?. In this domain the partition function is given by (see (5.25)) 



Z(J) 



QeGKi|i) 



dU exp 



Str 



z + J 
z 



iQ + Q-') 



, (6.22) 



where the integration is over the maximum super-Riemannian sub-manifold of 
GZ(1|1). This manifold is parametrized by 



Q = exp 



a 
P 







(6.23) 
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The integration measure is the Haar measure which in terms of this parameter- 
ization is given by (with 5Q = Q~^dQ) 



Sdet 



5Q 



kl 



5(j) Ss Sa SP 



dadpdcpds. 



(6.24) 



It is straightforward to calculate the Berezinian going from the variables 
{6Qn, to the variables {Sep, 6s, 6a, 6P}. The derivative 

matrix is given by 



B 



kl 



6(j) Ss Sa Sp 



1 I 

e^-''t' 
V e-'+'f ) 



2 


■s+i</> 



(6.25) 



Using the definition of the graded determinant one simply finds that SdetB = 
i. Up to a constant, the integration measure is thus given by d(f)dsdadp. 
We also need 



u _ ( cos(/.(l + ^f ) a{e' -e 



coshs(l — 



(6.26) 



After differentiating with respect to the source term {G{z) = dj log Z[J)\j=q) 
this results in (with x = VY^z) 



G{z) 



d(j)dsdadp a/? 



Cos,^(l + ^)eXcos^(l+if )-xcosh.(l-if )_ 

27r 2 



(6.27) 



With the Grassmann integral given by the coefficient of aP we obtain 



G{z) 



/ [cos + x{cos (j) + cosh s) cos (^Je^^'^"^ '^-'^^'^ . 



All integrals can be expressed in terms of modified Bessel functions. We find 

G{z) 



X 



Ii{x)Koix) + -il2ix)Koix) + Io{x)Koix) + 2Ii{x)Kiix)), 

(6.28) 



which can be further simphfied by the recursion relation hix) = Io{x) — 
2Ii{x)/x. As final result we obtain [51, 19, 20] 



G{z) 



x{Io{x)Ko{x) + h{x)Ki{x)). 



(6.29) 
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This calculation can be generalized to arbitrary Nf and arbitrary v. The 
calculation for arbitrary Nf with m j = is much more complicated, but with 
a natural generalization of the factorized parameterization, and using some 
known integrals over the unitary group, one arrives at the following expression 
in terms of modified Bessel functions 

^ = - + x{Ia{x)Ka{x) + Ia+l{x)Ka-l{x)), (6.30) 

V Zj X 

where a = Nf + This result is in complete agreement with the resolvent 
obtained [51] from integrating the microscopic spectral density (6.21). 

For a = this result is plotted in Fig. 3. We observe that, below some 
scale, lattice QCD data obtained by the Columbia group [58] closely follow 
this curve. The predictions of chRMT or of the partially quenched chiral La- 




Figure 3. 

The resolvent of quenched QCD. The points represent lattice data obtained by the Columbia 
group, and the theoretical prediction (6.29) is given by the solid curve. (Taken from ref. [45].) 



grangian have been studied by numerous lattice simulations [16, 17, 59, 50]. 
In all cases, agreement has been found in the expected domain of applicability. 
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7. Integrability and the QCD Partition Function 
7.1 Virasoro Constraints 

In this section we derive the small mass expansion of the QCD partition 
function by means of recursion relations for the partition function known as 
Virasoro constraints. The starting point is the QCD partition function in the 
ergodic regime given in (5.21). The quantities 

g,{tk) = det-' {M)Z,,Nf{M) (7.1) 

are invariant under the U{Nf) x U{Nf) transformations M ViMVf^. 
(Here and below we drop the superscript eff.) Therefore, Qu{tk) only depends 
on the eigenvalues of M'^ M which can be parameterized in terms of the mo- 
ments 

■ (7.2) 
A differential equation for the Qy{tk) is obtained from the unitarity relation 

Notice that the factor Y.V is included in M. This relation can be rewritten as 



1 d 



Qu{tk) = (7.4) 



_dM^adMlc dMc 

Using the chain rule and the assumption that the matrix elements of {MM^y~^ 
are independent for different values of s, we obtain [22] 

[Cs - Ss,i]g.{tk) = s>l. (7.5) 

The Virasoro operators Cg defined by 

= + W + (7.6) 

satisfy the Virasoro algebra 

[Cr, ^s] = {r - s)Cr+s- (7.7) 
Therefore, all Virasoro constraints with s > 1 are satisfied if Qi, satisfies 

[A-l]a^. = and C2g^. = 0. (7.8) 
This justifies the independence assumption above Eq. (7.5). 
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7.1.1 Solution of the Virasoro Constraints. We can expand Gi, as 

G^ = l + aih + 02*2 + autl + ■■■ . (7.9) 

From the first Virasoro constraint we obtain 

CiG„ = a2ti + {Nf + u){ai + 2auti) + ---, 

= 1 + aih + 02*2 + autl + ■■■ . (7.10) 

By equating the coefficients of the tj- we find 

ai = — ^— , a2 + 2(iV/ + i/)aii = ai. (7.11) 
i\ f -\- u 

From the second Virasoro constraint, C2G1, = aii+{Nf+u)a2 = 0, we obtain 
a2 = —aii/{Nf + u). Continuing this way we can obtain all coefficients in 
the expansion of the G{tk). This results in the small mass expansion of the 
partition function [60] 

ZQGD(M) ^ TVMMt 1 Tr(MMt)^ 
defM ^ ^ 4{Nf + u) ^ 32 {Nf + u){{Nf + u)^ - 1) ^ ^'^ ' 

An extension of this expansion to all three Dyson classes can be found in [61]. 

The small mass expansion can be used to obtain sum rules for the inverse 
eigenvalues of the Dirac operator [13, 60]. The QCD partition function can be 
expanded as (the prime indicates that 7^ 0) 

ZQCI^(M) = m'^^s(j^xf^{l + ^rf)) (7.13) 



k 



k 

This results in the expansion 

^u^^i^) = 1+^2^ /y-' QCD .7J4>) 

Imijn^Qm ""^^sZcr (M) '=A^ 
By equating this expansion to the expansion (7.12) for equal masses given by 

4(iV/ + v) 

we obtain a Leutwyler-Smilga sum rule [13] for the inverse Dirac eigenvalues 
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7.1.2 Flavor-Topology Duality. If we construct an {Nf + z/) x {Nf + u) 
matrix M with Mij = Mij for z, j < Nf and Mij = otherwise, we have that 
Tr (MM+)^ = Tr(MMt)fc. Since the Virasoro constraints only depend on the 
combination Nf + v we have 

det-^MZ.^NfiM) = Z,=o,Nf+u{M). (7.17) 

This relation is known at the flavor-topology duahty [62]. 

7.2 T-Function 

The unitary integral in the QCD partition function can actually be evaluated 
analytically for an arbitrary number of flavors. We will show that it can be 
rewritten as a r-function. The unitary matrix integrals can then be evaluated 
by means of a Harish-Chandra-Itzykson-Zuber type integral and the use of 
flavor-topology duahty. 

7.2.1 Itzykson-Zuber Integral. We consider the integral 

1 = j dC/dFe5^(^'«^^+^^'^^), (7.18) 

where U G U{Ni) and V G U{N2)/U^^{1) and the integral is over the Haar 
measure of these groups. The matrices R and S are arbitrary rectangular com- 
plex matrices. Without loss of generahty, they can be taken diagonal with 
Rkk = rk > and Skk = Sk > and all other matrix element equal to zero. 
Using the diffusion equation method one can derive the result [63] 

A{{sl})A{{rl}y ^'-'^^ 

This result first appeared in the Russian hterature [64] as a solution of the 
Laplace equation. It was proved independently in [63]. 

7.2.2 The QCD Partition Function is a r-Function. In this subsection 
we show that the finite volume QCD partition function is a r-function. Using 
the flavor-topology duality (7.17) with M^k = Xk for k < Nf we can write 

Z,iv,(M) = de™ / ^Ue^MMU^+MW). (7.20) 

' JueU{Nf+u) 

This integral is exactly (7.18) with R equal to the Nf x [Nf + u) matrix with 
rk = Xk- A finite result is obtained if the Nf diagonal matrix elements of S 
are expanded as = 1 -I- Ssk- For Ssk Ov/e obtain 

I.i^ksi) = ^-^/^-^)(xfe)(5sO^-^^^^fc,S,7. (7.21) 
.7=1 .1=1 
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with (the upper index between brackets such as (fc) denotes the /c'th derivative) 
Akj = 4-'4^-'Hxk), Bji=^^0-^. (7.22) 

Up to a constant, the determinant of B is given by 

detB ^ A{{{l + 5skf}), (7.23) 
and cancels against the denominator in (7.19). We finally obtain the result 

Z.,N,{M) - . (7.24) 

This results was first obtained in [65] and independently for equal masses in 
[9, 66]. Using identities such as 

x^dl = {xd^f - xd^, x^dl = {xd^f - ?,{d^f + 2xd^, (7.25) 

we can rewrite this partition function in terms of derivatives 5^ = x^dx^ as 

^^,7V/ (M) = ^(i^U^ ■ (7-26) 

This form of the partition function is also known as a r-function [22, 67]. 

7.2.3 QCD Partition Function for Equal Masses. The limit of equal 
masses in the partition function (7.26) can be obtained by writing 

Xk = x{l + 5xk), (7.27) 

and taking the limit 8x}^ — ^ 0. Because all columns are the same for (5xfe = 0, 
we have to expand the matrix elements to order {5xk)^f~^ . The expansion to 
this order can be combined into 

del[5i-h,{xk)] = dei[x'-\5i-^h{x)f-^^] det[^^^^]. (7.28) 

Using that 

det[(fe,)'-^] = A({fefe}), A({4})=a;^/(^/-i)A({fefc}),(7.29) 

we obtain the partition function 

Z^^Nf{x) = cx-^f^^f-^^ det[x^-\{xdxy-^I,.{x)f-% (7.30) 

With the help of the identities (7.25) the derivatives can be combined into 
derivatives {xdxY resulting in 

Z,,Nf{x) = cx-^/(^/-i) det[(xa^)'+^-27,(x)]i<,-Kiv,. (7.31) 
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7.3 Toda Lattice Equation 

In this section we show that the QCD partition function for equal masses sat- 
isfies a Toda lattice equation. This result and generalizations thereof were first 
obtained in [68]. The Toda lattice was originally introduced as a one dimen- 
sional lattice in which neighboring atoms interact via a potential that depends 
exponentially on the distance. The Hamiltonian equations of motion of this 
system can be written in the form of the Toda lattice equation discussed below. 
Because of the existence of a Lax pair, they have infinitely many constants of 
motion. For a more elaborate discussion of the Toda lattice equation and the 
relation to integrable systems, we refer to [29, 69, 70]. Several subsections 
below are based on the paper by Forrester and Witte [71]. 

7.3.1 The Sylvester Identity. We all know how to expand a the deter- 
minant matrix with respect to its co-factors given by 

Cij = ^T-deiA. (7.32) 

What is less known is that there exists a remarkable identity that relates co- 
factors to the double co-factors defined by 

This identity, which is known as the Sylvester identity [72], is given by 

CijCpq CiqCpj = det ACij-pq. (7.34) 

For example, it holds for a 2 x 2 matrix with i = j = 1 and p = q = 2. 

7.3.2 Toda Lattice Equation. We apply the Sylvester identity to the 
determinant that appears in the partition function (7.31). For i = j = Nf — 1 
and p = q = Nf v/e obtain 

CNf-l,Nf-lCNf,Nf — CNf-l,NfCNf,Nf-l = det ACNf-l,Nf-l;Nf, Nf , ■'^5) 

with matrix A given by 

Aji = {xdJ+^-^I,{x). (7.36) 

The derivative of a determinant is equal to the sum of determinants with one of 
the rows replaced by its derivatives, or is equal to sum of the determinants with 
one of the columns replaced by its derivative (in both cases we have in total Nf 
terms). For the matrix A in (7.36) only differentiating the last row or column 
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gives a nonzero result. This allows us to rewrite the co-factors as derivatives 
of det A. In particular, we find 

CNf-i,Nf = —xd^det A]sr^_i, Cn^^n^^i = —xd^det A]sf^_i, 
CNf,Nf = det^iVj.-!, CNf-i,Nf-i = ixdx)'^detANf-i- 

CNf-l,Nf-l;Nf,Nf = detANf-2- (7.37) 

To obtain the second last identity we first differentiate the columns and then 
the rows. Inserting this in (7.35) we find 

dct^Tv, det^Arf-2 

(x5^)2logdet^iv,-i = T^^-^ (7.38) 

det An^^i 

Next we substitute the relation between det Anj, and the partition function 

det^AT, = -x^/(^/-i)zf^(x). (7.39) 

^ c 

The prefactor contributes a factor to the r.h.s. of (7.38) but does not con- 
tribute to its l.h.s.. After raising Nf by 1 we obtain the celebrated Toda lattice 
equation [22, 67] 

{xd.flogZ.M-) = cx-^-El±p^p^^. a AO) 

['^i',Nf{X)\ 

The large-x limit of the partition function for equal masses, 

Z.N fix) = I dU dot'' Ue^^^^+^~'\ (7.41) 

JueU{Nf) 

is obtained by performing the [/-integral by a saddle point approximation in- 
cluding the Gaussian fluctuations. The asymptotic behavior is given by 

^Nfx 

Z.,Nf{-)^^^. (7.42) 
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Using this result to normalize the partition function we find that c = Nf. 

7.4 Painleve System 

In this subsection we derive the Toda lattice from the Backlund transforma- 
tion of a Painleve equation [71]. 

The partition function (7.41) can be obtained from the "double" scaling limit 
of the random matrix partition function 

k=l 

(7.43) 
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where x = niNT, is kept fixed for N ^ oo. Using = A|S^A^^ + as 
new integration variables, we obtain in this limit 



^™(-) = T{[rdxl{xl-xyxf^V^^e^''l^e-T..^l' 

7- JX^ 



1/2N 

k -""^ 

= x''^fe''^/'^EN{[0,s = x^],Nf,i^). (7.44) 

En{[0, s = x'^],Nf, v) can be interpreted as the probabiUty that there are no 
eigenvalues in the interval [0, .x^] for the joint probability distribution given 
by the integrand of (7.44). For v = this is the probability of the partition 
function with topological charge Nf and no flavors. If we introduce a{t) by 

E!,{[0, s], Nf, z.) = e-/o f (-(*)+iKiV/+-))^ (7.45) 

then the function a{t) satisfies the Painleve equation [71], 

ita"f - {N} - iy^){a'f + a' {4a' - l){a - ta') - ^ = 0. (7.46) 

The boundary conditions for this differential equation follow from the asymp- 
totic behavior of the partition function (7.41). Using (7.42) we find for the 
large x behavior of En{[0, s = x^,Nf, u), 

Er,{[0, s = x\ Nf, v) ~ ^-^/V2--iV/giV,.-.V4^ (7_47) 

so that the large-s behavior of a{s) is given by 

a{s)^^^s + --- + ^. (7.48) 
The Painleve equation (7.46) can be derived from the Hamiltonian [71] 

tH = qY - {q^ + {Nf + iy)q - t)p + iV/g (7.49) 
with the identification 

a{t) = - {tH)\^^, - ^{Nf + u)u + i (7.50) 

For example, we have the Hamiltonian equations of motion 

{my = p, {tH)"=p' = -^, tq' = ^. (7.51) 

Such Hamiltonians play an important role in the theory of exactly solvable 
models. Hamiltonians with different values of u and Nf are connected by 
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a Backlund transformation. This is a canonical transformation together with 
{u, Nf) {i>, Nf) such that, in the new variables, the same Painleve equation 
is satisfied. In our case, we have the Backlund transformation (at fixed 

T : Nf ^ Nf + 1, 

HNf Hn^+i = Hnj: + qNf — QNfPNf, 

QNf QNf+l, 

PNf PNf+1- (7.52) 

Below we do not need the explicit transformation rules of gjv^ and pN^, but of 
the inverse transformation of qNf {pNf — 1) which is given by 

T-^Nf ipNf - 1) ^ -qNf {pNf-l)+Nf + v-^. (7.53) 

PNf 

If we define the r-function by 

TNf=J'^^''f'^\ (7.54) 
one can easily derive the equaUties 

tdt log — - — 2 — = tH]sff-i — tH^f + tH^f+i — tH^f, 

'^Nf 

= -T~^{qNf{'i--PNf) + {qNf{'i--PNf), 

= tdtlogdt[tHNf], 

= tdt log dt [tdt log TNf ]. (7.55) 

To derive the second last equality we have used the inverse Backlund transfor- 
mation and the Hamilton equations (7.51). Integrating this equation once and 
putting the integration constant equal to zero we find the Toda lattice equation 

{tdtflogrNf=t^-^^^i^^. (7.56) 

7,4.1 Solutions of the Painleve equation. The probability E{[0, s],Nf, u) 
is related to the partition function (7.41) by 

E{%s\,Nf,u) = s-'^f'^/^e-'/^Z^^NfiV^). (7.57) 
For Nf = the partition function is normahzed to 1 so that 

2 
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Indeed this is a solution of the PI 1 1 Painleve equation (7.46). 
For Nf = Iv/e have that 

E{[0, s], Nf = l,u) = s-^/2e-^/^I,(Vi) = e-/o f (7.59) 
resulting in another solution of (7.46): 

s iP' d 

crNf=i{s) = 4 ~ y ~ ^og I^{^/s). (7.60) 

For a; — > the modified Bessel function behaves as Ii,{x) ^ so that the 
t-integral in (7.59) is well-behaved for t ^ 0. Only recursion relations for 
Bessel functions are required to show that (7.59) is a solution of the Painleve 
equation. Since {—lYK^{x) satisfies the same recursion relations as Iu{x) 
this provides us with another solution of the Painleve equation. This solution 
corresponds to the partition function with Nf = —1 where 

E{[0,s],Nf = = s-'/^e-'/'K^iVs) = e'/o f (-W+|-(-i+-))^ (7.6I) 

and satisfies the boundary condition for Nf = —1. 

7.4.2 The Bosonic Partition Function. The natural interpretation of 

Nf = — 1 is as a bosonic flavor. In this section we will derive the low en- 
ergy limit of the QCD partition function for Nf bosonic flavors with equal 
masses. We already have seen in section 3.3.2 that the Goldstone manifold for 
n bosonic quarks is given by Gl{n)/U{n). Using the same invariance argu- 
ments as before one obtains the low-energy effective partition function 

Z,_„= / det'^(Q)e^^^^^(«+«"'). (7.62) 

jQeGl{n)/U{n) 

In this case Q can be diagonalized as Q = [7diag(e^'=)[/^^, so that an eigen- 
value representation of this partition function is given by [73] 

/ n^^^^n^'''' n(e'' - e'0(e""= - e-")e''^k^'^'K (7.63) 

k k k<l 

The Vandermonde determinant can be written as 

[Ke^/c - e^O = det[eP'^]o<p<n-i,i<q<n (7.64) 

k<l 

and a similar expression for Sk — Sfc. By expanding the two determinants the 
integrals can be written as modified Bessel functions which can be combined 
into a determinant as follows [73] 



Z„-n{x) = C-ndet[Ki,+k+i{x)]o<k,l<n-l- 



(7.65) 
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From the observation that {—iyK^{x) and Iu{x) satisfy the same recursion 
relations, and that the factor (— 1)*^ does not affect the determinant, (7.65) can 
be rewritten as the r-function 

Z,,_nix) = det[(xa^)'=+iZ,,_i(x)]o<fc,Kn-i. (7.66) 

with Zy^-\(x) = Ki,{x). The bosonic partition function can also be analyzed 
along the same Unes as the fermionic partition function. On the other hand 
this derivation can be simply modified to obtain the partition function for Nf 
fermionic flavors with equal mass. 

The bosonic partition functions thus satisfy the same Toda lattice equations 
as the fermionic partition functions. The semi-infinite hierarchies are con- 
nected by 

lim -{xd^f log Z^^ri{x), (7.67) 
which is related to a derivative of the resolvent. 

7.5 Replica Limit of the Toda Lattice Equation 

The resolvent can be obtained from the repUca limit 

G(z) = lim ilogZ^,„(z). (7.68) 

n^O n 

If we take the replica limit of the fermionic (n > 0) or bosonic (n < 0) 
partition functions directly, we will obtain a result that differs from the super- 
symmetric calculation. These problems can be avoided if the take the replica 
limit of the Toda lattice equation (7.40). With the normaUzation Zi,^q{x) = 1 
we obtain the relation 

xd^xGix) = 2x^Z^,i{x)Z^,_i{x). (7.69) 

Inserting the expressions for Zi,^i and ^1^,-1 we find [23] 

G{x) = -+ x{K^{x)U{x) + (7.70) 
x 

which agrees with the result obtained from the supersymmetric method (6.30). 
This result has also been derived from the solution of the Painleve equation 
(7.46) forn ^ [24]. 

The validity of the replica limit of the Toda lattice equation can be proved by 
extending the Toda lattice hierarchy to include an additional spectator boson 
with mass y and using the identity [74] 

\im -{xdxf log Ziy^nix) = lim xdx{xdx + ydy) log y). (7.71) 

n — 77/ y 
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7.6 Replica Limit for the GUE Two-Point Function 

We have two possibilities for the generating function of the two-point func- 
tion of the Gaussian Unitary Ensemble: a fermionic generating function or a 
bosonic generating function. The fermionic (bosonic) generating function for 
the two-point function is defined by 

Zn{x, y)= J dHP{H)det''{x + ie + /i")det"(y -ie + H), (7.72) 

with n > (n < 0). We will consider the microscopic limit where ttIx — 
y)Np[x) = r is kept fixed for ^ oo. In that case the two-point function in 
the center of the spectrum only depends on r and is given by 

R2{r) = - lim ^djZ,,{r), (1.13) 

n— >0 n 

both in the fermionic and the bosonic case. In an eigenvalue representation 
of the Goldstone fields the microscopic Umit of the generating function Zn{r) 
can be written as [21, 37] 

^n{r) = C ^duk \[{uk - uife'^^^-K (7.74) 
k k<l 

This partition (7.74) can be written as a r-function. The first step is to expand 
the Vandermonde determinant 

Next we use that 

J^^ dukule''-''^ = {dirfZ^ (r) , (7.76) 

which results in [75] 

Z„(r) =n![det(5i,)^+^Zi(x)]o<i,,<n-i. (7.77) 
The partition function Z\ {x) is given by 

Zi(r) = due'''''. (7.78) 

The microscopic limit of bosonic partition function can be rewritten sim- 
ilarly. The main difference is the convergence requirements of the bosonic 
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integrals which are essential for the structure of the Goldstone manifold. In an 
eigenvalue representation of the Goldstone fields we find [21] 

poo „ 

Z_„(r) = / Y[duk Y[{uk - uife'^-iZk^K (7.79) 
■'^ k k<l 

This partition function can also be written as a r-function. By expanding the 
Vandermonde determinant we can express this generating function as a deter- 
minant of derivatives 

Z_n{r) = n\[detidiry+^ Z_i]o<i,j<n-i, (7.80) 
with Z_i (r) given by 

Z_i{r)= / due'""^. (7.81) 



dir log Zn (x) = n , (7.82) 



1 

Because of the derivative structure of the partition function, we can again 
use the Sylvester identity to derive a Toda lattice equation. In this case we find 

^2 i„„ ^ _ „2^n +l(^)^«- l(0 

where the factor n? follows from the choice of the normalization constants. 
We have made this choice because the left hand side is proportional to n^. The 
two-point correlation is given by the replica Umit of (7.82) 

R2{r) = -lim^a2logZ„(r) = Zi(r)Z_i(r) 

= t due''''' r e'"^ = 2i^-^ — , (7.83) 

J-l Jl XX 

which is the correct analytical result for the two-point function. This deriva- 
tion explains the factorization of the two-point function into a compact and a 
non-compact integral which characterizes the result obtained by a supersym- 
metric calculation [76]. The fermionic partition functions, the bosonic partition 
functions and the super-symmetric partition function form a single integrable 
hierarchy which are related by the Toda lattice equation [23]. A closely related 
way to derive the two-point function of the GUE is to take the replica limit 
of the corresponding Painleve equation. For a discussion of this approach we 
refer to [24] which preceded our work [23] on the Toda lattice. 

8. QCD at Finite Baryon Density 

In this Chapter we study the quenched microscopic spectrum of the QCD 
Dirac operator at nonzero chemical potential when the Dirac operator is non- 
Hermitian with eigenvalues scattered in the complex plane. Using the replica 
limit of the Toda lattice equation we obtain the exact analytical result for the 
microscopic spectral density [25]. 
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8.1 General Remarks 

The average spectral density of a non-Hermitian operator is given by 

p{X) = {J2s\X-Xk)), (8.1) 

k 

and the average resolvent is defined by 

Using that dz*{l/z) = Tr6'^{z) we easily derive 

d,*G{z)U=x = MX). (8.3) 

The resolvent can be interpreted as the electric field in the plane at point z from 
charges located at the position of the eigenvalues. For example, Gauss law is 
given by 

^ G{z)dz = 2mQ, (8.4) 
where Q is the number of eigenvalues enclosed by C. 

8.2 The Ginibre Ensemble 

The obtain a better understanding of the resolvent for a non-Hermitian ran- 
dom matrix ensemble, we first consider the Ginibre ensemble [77] defined by 
the probabiUty distribution 

^(C) = e-^^^^\ (8.5) 

with C a complex N x N matrix. The eigenvalues of C wee given by the solu- 
tions of the secular equation det(C — A^) = 0. If all eigenvalues are different, 
the matrix C can be decomposed as 

C = VAV-^, (8.6) 

where F is a similarity transformation and A = diag(Ai, • • • , Aat). The joint 
eigenvalue distribution is obtained by using the decomposition 

C = UTU-^, (8.7) 

with U a unitary matrix and T an triangular matrix with the eigenvalues of C 
on the diagonal. After integrating out the upper triangular matrix elements we 
obtain (see the lectures by Zabrodin [78] for a derivation), 

p({Afc}) = |A({Afc})pe-^Ej^fc^ (8.8) 
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This distribution can be interpreted as repulsive charges in the plane balanced 
by an external force N\z\. The resolvent is equal to the electric field of the 
eigenvalues. For an equilibrium distribution we have G{z) = N\z\. Using that 
the eigenvalue density is spherically symmetric, we find from Gauss law 

27rr|G| = 27r r p{r')dr', (8.9) 
Jo 

so that p{r) = N/tt. Because the total number of eigenvalues is equal to N, 
they are located inside the circle \z\ = 1. The resolvent is thus given by 

TV 

G{z) = Nz*e{l-\z\) + —e{\z\-l). (8.10) 
8.3 QCD at Nonzero Chemical Potential 

The QCD partition function at nonzero chemical potential n is given by 

ZQCD=Ee"^^'''"'^'^'=^ (8.11) 
k 

where Ef^ is the energy of the state, and is the quark number of the state. 
At zero temperature (/3 oo) the partition function does not depend on for 
/X < m]sf/N]sf, where N is the particle with the smallest value of rriN/Nj^j. 
For QCD A'^ is the nucleon with quark number = 3. This implies that the 
chiral condensate does not depend on for /x < uin/Nn. 

The QCD partition function can be written as a Euclidean path integral with 
the fermionic part of the Lagrangian density given by 

jC = -ipDip + m'ip'ip + fj,'tp^oip. (8.12) 

with D the anti-Hermitian Dirac operator Since //70 is Hermitian, the Dirac 
operator as a whole is non-Hermitian. As a consequence, the eigenvalues are 
scattered in the complex plane [79]. The fermion determinant is in general 
complex. This means that it is not possible to study the QCD partition function 
by means stochastic methods which severely Umits our knowledge of QCD at 
nonzero chemical potential. 

The question we wish to address is if there is a domain where the fluctu- 
ations of the Dirac eigenvalues are universal and can be obtained from a ran- 
dom matrix partition function with the global symmetries QCD, or equivalently 
from a chiral Lagrangian. In this domain we will calculate the resolvent and 
the spectral density from the replica limit of the Toda lattice equation [25, 26]. 

8.3.1 Generating Function for the Quenched Spectral Density. The 

quenched spectral density is given by the replica limit [80, 81, 82] 

pi"^°(z, z*,fi) = lim —d,d,* log Z,,niz, z*,p), (8.13) 
n— >o 7rn 
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with generating function given by 

z\ii) = (det"p + /no + z)det"(-Z) + ^i^^ + z*)). (8.14) 

The product of the determinants in (8.14) can be written as the determinant of 

[83, 84] 



I id + IX z 0\ 

id-n z* 

z id^ + n 

V z* id^ -n J 



id + /xi MflL 
iS + H2 



(8.15) 



Here we have used the decomposition of the Dirac operator given in (3.9). We 
observe that the U{2n) x U{2n) flavor symmetry is broken by the chemical 
potential term and the mass term. Invariance is recovered by transforming the 
mass term as in the case of zero chemical potential (see (5.2)) and the chemical 
potential term by a local gauge transformation [85]. For the chemical potential 
matrices the latter transformation is simply given by 

W ^ URfiiU^\ ^l2 ^ ULfJi2UL^. (8.16) 

The low-energy limit of quenched QCD should have the same transformation 
properties. In the domain /x <C l/L and zS <C F"^ jl?, we only have to 
consider the zero momentum modes. Using that the Goldstone fields transform 
as [/ — > UrUU^^ we can write down the following invariants to first order in 
the quark mass and to second order in the chemical potential 

Tvfil TTU-^fiiUfi2, TxMrlU-\ TtMlrU. (8.17) 

The low energy effective partition function is therefore given by 

Z,r^{z,z*,^^) = I dC/det'^C/e-^^[^'«l[^"'^l+^^^(^+^"),(8.18) 

Ju{2n) 



5=( n" ^ V M=(']: ° ). (8.19) 



where 



-In J' V Z*ln 

8.3.2 Random Matrix Model. The partition function (8.18) can be 
obtained from the large limit of a random matrix model with the global 
symmetries of the QCD partition function. For = 0, the model is defined by 
an integral over iV/2 x complex matrices [14, 86, 82], 



Z^=o{mf,fi) = dWjl det(£'(Ai) +m/)e-^^"^^^\ (8.20) 
J f=i 
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The Dirac matrix has the structure 

iW + i2 



For QCD with three or more colors in the fundamental representation, the ma- 
trix W is complex ((3d = 2). One can also introduce random matrix ensembles 
with /3£) = lor/3£) = 4by choosing the matrix elements of W real or quater- 
nion real, respectively [87]. 

An alternative random matrix model [27] is obtained by replacing the iden- 
tity matrix that multiplies by a complex matrix with the same distribution as 
W. This random matrix model is in the same universality class but turns out to 
be mathematically simpler In particular, the joint eigenvalue distribution has 
been derived [27] which makes it possible to calculate correlation functions by 
means of the orthogonal polynomial method. 

8.3.3 Mean Field Analysis. The macroscopic spectral density of the 

partition function (8.18) can be easily obtained by means of a saddle point 
approximation [84]. Using an Ansatz that is diagonal in replica space, 

( cos 9 e*<^sin6' \ 
^-[-e-i't>sme cose J' ^^"^^^ 

the partition function is given by 

The extrema are at 

cos0 = l, or cos0=^^y^. (8.24) 
The critical value of /i is at the point where the two saddle points coincide, 

= (8.25) 
The partition function at the saddle point is given by 

/x>/.e : Z„(z,z*,M) = e-^(^'^'^'+^'(^+^*)V(«^'^')). (8.26) 
For the resolvent and the spectral density we thus find 

^i<^ic : G^^''\z,z*,^,) = V^, p'^^''''{z,z*,fi) =0, (8.27) 

vy-^(7 -I- 7*) V'^v 



The eigenvalues are located inside a strip of width 4F^/x^/E. 
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8.4 The Microscopic Spectral Density 

The microscopic spectral density cannot be obtained from a mean field anal- 
ysis. The assumption in the derivation of the previous section is that saddle 
point is proportional to the identity in replica space, so that the replica limit 
can be obtained from the calculation with one rephca. The generating function 
for the microscopic spectral density depends in a nontrivial way on the number 
of replicas which, as we have seen before, can be obtained from the replica 
limit of a Toda lattice equation. In this subsection we closely follow [25]. The 
starting point is a remarkable integration formula to be discussed next. 

8.4.1 Integration Formula. By decomposing a U (2n) matrix as 

ui \ f vi \ f Vl — ft^ b 1 / '^1 



»Jl .J[ 6 -VT^)[ ' .tj-»^« 

with ui, U2, Vi G U{n), V2 G U{n)/U"'{l) and 6 a diagonal matrix, the fol- 
lowing integration formula can be proved [25] 

JU{2n) 

where c„ is an n-dependent constant and 

Z,,,{x,y) = [' AdA74Ax)J.(-Ay)e'S.«^^°^('^^°^"''\ (8.30) 
Jo 

8.4.2 Toda Lattice Equation at Nonzero Chemical Potential. Using 
the integration formula (8.29) for p = 1 we find that the zero momentum 
partition function Zi,^n{z, z*,^) (see Eq. (8.18)) can be written as 

Z,,n{z,z\ii) = det[(za,)^(z*a,.)'^.,i(^,^*)]o<fc,Kn-i(8.31) 
where 

Z^i{z,z\^j) = /^AdAe-2^^''^'(^'-i)|I^(AzyE)|2. (8.32) 
Jo 

By applying the Sylvester identity to the determinant in (8.31) for i = j = 
n — 1 and p = q = n and expressing the cof actors as derivatives, we find a 
recursion relation that can be written in the form of the Toda lattice equation 

o *o , 7 / * N 2 Zu,n+iiz, z*,n)Z^^ri-i{z, z*,ii) 
zd^z d^^ log Z^^n{z,z ,fj,) = —{zz ) [Z^^{z z* n)]^ ^ ^ 
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For the spectral density we find tlie simple expression {Zi,fi{z, z*,iJ,) = 1) 

1 zz* 
p^'^^^z.z^iJi) = lim—dAAogZ,,r.iz,z*) = —Z,,i{z,z*)Z,,_i{z,z*). 
n— >o Trn z 



(8.34) 



What remains to be done is to calculate the bosonic partition function for n = 
— 1 which will be completed in the next subsections. 



8.4.3 The Bosonic Partition Function. In this subsection we evaluate 
the low-energy limit of the QCD partition function at nonzero chemical poten- 
tial for one bosonic quark and one conjugate bosonic quark and u = 0. We 
closely follow [25]. Because of convergence requirements, the inverse deter- 
minants of nonhermitian operators have to be regulated. This is achieved by 
expressing them as the determinant of a larger Hermitian operator [83] 



det" 



z id + 11 
iS + z 

( 



det 



limdet ^ 

e-»0 



/ 



N/2 



_i / z* —id + fi 

\ -id^ + fi z* 

e z id + fi \ 

e id^ + 1^ z 

z* —id + jjL e 

V -zrft + 11 z* e y 

/ e z id + /J, 

z* e id — iJ, 

-id^ + ljL e -z* 



V 







(8.35) 



4] 



The mass matrices are given by 



Ci = 



€ Z 

z e 



and 



C2 



fj, —z 



-z^ 
e 



kl 



-ICiI. (8.36) 



with / = i(T2. For the random matrix model (8.21) we have that d = W. 
The Gaussian integral over W results in the 4-boson term exp[— 2Tr(5i Q2/N] 
with 



^1 • 92 
^*2-<t>2 



Q2 = 



^3 ■ </^3 (P3- </^4 



(8.37) 



Qi = 

and we have used the notation 
Hubbard-Stratonovitch transformation, we linearize the 4-boson interaction 
term by the Hermitian matrix 5 function 



'^k ' '^i = ^i^=i ^k*^l- Instead of the usual 



(8.38) 
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where the integral is over Hermitian matrices F. We thus find the identity 



exp 



-;|lVQi Q2 



(8. 



The integral over the cpk is uniformly convergent in F and G which justifies the 
interchange of the order of the integrals. This results in the partition function 

Z_i = ldQidQ2ldFdGe^'^-''^^'^'-''^''^^+''^'^^^^'-'^''^-^'^''^'^ 

xdet-fr^ + ^ ^-l,), (8.40) 
\ e + G J ' 

where we have used a block notation and the mass matrices (8.36). We have 
also simplified this integral by changing integration variables according to 
F ^ F-Cl andG ^ G + ICfl and Qi NQi/2, i = 1,2. For reasons of 
convergence we have kept the infinitesimal increments inside the determinant. 
In the weak nonhermiticity limit, where fi^N is kept fixed for N — > 00, the 
determinant can be approximated by 

'^^*~^('/if J+g) " det-T(e + F)det-T(e + G) (8.41) 

(1 



X exp 



Nfi^ 1 1 
^^^^^^^ 



The F and G variables in the jj.'^N term can be replaced the saddle point values 
of F and G at = given by {e+F)Qi = i and {e+G)Q2 = i. The remaining 
integrals over F and G are Ingham-Siegel integrals given by [88] 

j dFdet-'^(e + F)e!-^'^^ - e(Q)det"-f (Q)e-^^^'5, (8.42) 

where the integral is over pxp Hermitian matrices, Ime < 0, and 9{Q) denotes 
that Q is positive definite. These manipulations result in 

Z_i{z,z\ii) = J dQidQ2e{Qi)9iQ2)ded-\QiQ2) (8.43) 



xe 



T^[ifCliQl-IQ2l)-fQlQ2-^t^^Qlcr3Q2Cr3] 



In the limit AT 00 the integrals over the massive modes can be performed 
by a saddle point approximation. The saddle point equations are given by 

Qr^ -Q2 = 0, Q2 ^ - Qi = 0. (8.44) 
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Both equations can be rewritten as 

Qi = Q2\ (8-45) 

and therefore only four of the modes, which we choose to be Q2, can be in- 
tegrated out by a saddle-point approximation. The quadratic fluctuations give 
rise to a factor 7r^/det^Qi. The integral over the remaining modes has to be 
performed exactly. We thus arrive at the partition function [25] 

Z_i(z,z*,m)= / ^^0(Qi)e^[^fCf{Qi-^Qr^J)-fM^Qi-3QrV3].(8.46) 
J det Qi 

Before evaluating this integral, we rederive this partition function based on the 
symmetries of the QCD partition function. 

8.4.4 Symmetries of For jj. = and Ci = C2 = the sym- 

metry of the partition function (8.35) are the Gl{2) x Gl{2) transformations, 



2 / V '^2 / \ 92 / \ 92 



01 \ 


f 0t 


02 J ' 


I 02 


03 A 


f 0,3 


04 ; ' 


I 04 



where we have disregarded convergence. This symmetry can be extended to 
nonzero mass or chemical potential if we adopt the transformation rules 

Ci ^ U2(iU{\ C2 ^ UiC2U2\ 

/ii ^ U2iiiU^\ M2 ^ UiH2U{\ (8.48) 

where //i is the chemical potential matrix that is added id and /Lt2 is the chem- 
ical potential matrix that is added to —id). These matrices are introduced for 
the sake of discussing the transformation properties of the partition function 
(8.35) and will ultimately be replaced by their original values //i = At2 = l^o^. 
The chiral symmetry is broken spontaneously to Gl{2) by the chiral conden- 
sate. Because the bosonic integral has to converge, the Goldstone manifold is 
not Gl{2) but rather Gl(2)/U(2), i.e. the coset of positive definite matrices as 
in the case of zero chemical potential. Under a Gl{2) x Gl{2) transformation 
the Goldstone fields transform as 

Q ^ UiQU^\ (8.49) 

The low energy effective partition function should have the same transforma- 
tion properties as the microscopic partition function (8.35). To second order in 
H and first order in the mass matrix we can write down the following invariants 

TVCiQ, TVC2Q~\ TrQ/xiQ-V2, Trfij, Trul (8.50) 
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We also have the discrete symmetry that the partition function is invariant un- 
der the interchange of Ci and ^2- This symmetry implies that the coefficients 
of the two mass terms in the effective partition function are the same. Using 
that the integration measure on positive definite Hermitian matrices is given by 
dQ/det^Q, we finally arrive at the effective partition function 

Z,Uz,z*)= I A^^''{Q)dQ ^_EhpLT.\QM\Q-\B]+^^Tr(C,Q+C2Q-^) ^ 

jQeGi{2)/u{2) det^{Q) 

(8.51) 

The partition function (8.46) for = is recovered after making the identifi- 
cation y ^ AT, E ^ 1 and F2 ^ 1 and Q ^ Q^. 

8.4.5 Calculation of the Integral over Q. To evaluate the integral 
(8.51) we use the parameterization 

p. _ t ( ^ cosh s e*^ sinh s \ ,„ 
I e~'^sinhs e~''coshs j ' 

where 

rG( — oo,cx)), sG (—00,00), te(— 00,00), G (0, vr) . (8.53) 

The Jacobian relating the measures dQ/det^Q and drdsdtdO is given by 

J = 4e^* cosh s sinh s . (8.54) 

We first perform the integral over r, which gives a factor 2Ko {2Ne cosh s cosh t) 
with leading singularity given by ~ — log e. This factor is absorbed in the nor- 
mahzation of the partition function. Then the integral over 9 gives a Bessel 
function. Introducing u = sinh s as new integration variable we find [25] 

/oo roo 
dt / due^^^Joi^Vuix^ cosh^ t + sinh^ tf^) 
-00 Jo 

To do the integral over u we use the known integral 



da;a;"+^e-"^V„(/3x) = ,^ ^° . e"^'/^". (8.56) 



lo ^"^"^ ^ ""'"'■^^ (2a)»+i 



This results in 



n ^p-Vfi^F^ poo _^ ^ 



(8.57) 
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Using that cosh^ ^ = 5 + I cosh 2t and sinh^ ^ = ~ 5 + I "^o^^ 2t, the integral 
over t can be rewritten as a modified Bessel function resulting in [25] 

= x + iy, z*) = C_i e 4,^-^ [ 1 • (8.58) 

8.4.6 The Dirac Spectrum at Nonzero Chemical Potential. The final 
result for the quenched spectral density is obtained by substituting the partition 
functions -^{z, z* , fi) and Z^^_i{z = x + iy,z*,^) in expression (8.34) 
obtained from the replica limit of the Toda lattice equation. We find, 



X /'AdAe-2^^'A^'^')|/,(AzyS)|2. (8.59) 
Jo 

The normalization constant has been chosen such that the n ^ limit of 
^qucn^^;^ y, fi) for large y is given by T,V/tt (see below). 

In the limit Re{z)T, / jj? F"^ <C 1 the upper limit of the integral in (8.59) can 
be extended to infinity. Using the known integral 

the spectral density can be expressed as 

2 VTi"^ 
p'i'''"'{x,y,n) = -u^zz*K^{zz*u)I^{zz*u) with u = —r-.. (8.61) 

Therefore, the spectral density becomes a universal function that only depends 
on a single parameter u. This parameter can be rewritten in a more physical 
way as u = 7rp'^'*^™(a;, y, ji). For the dimensionless ratio we obtain 

„quen(' \ 

f- yjhrL = 2uzz*K^izz*u)U{zz*u), (8.62) 

which is universal combination that depends only on a single universal combi- 
nation zz*u. (This result was obtained in collaboration with Tilo Wettig). 

In the thermodynamic limit the Bessel functions can be approximated by 
their asymptotic limit. This results in 

= / dXe~^^^ . (8.63) 
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For y — >^ oo the integral over A can be performed by a saddle point approxi- 
mation. If the saddle point is outside the range [0, 1] the integral vanishes for 
V ^ oo. We thus find for the spectral density 

p^"-(x,y,/.) = ^-^ for (8.64) 

and p^^'^'^{x,y, ji) = outside this strip. This result is in agreement with 
the mean field analysis [84] of the effective partition function given in section 
8.3.3. For the integrated eigenvalue density we find 



cZxpl"^°(ar, y,n) = — (8.65) 



in agreement with the eigenvalue density at = 0. 

9. Full QCD at Nonzero Chemical Potential 

In quenched QCD the chiral condensate G(m, fi) m/jj? in the region 
where the eigenvalues are located. In section 8.3 we have argued that that 
chiral condensate in full QCD does not depend on fi for jj, < itin/Nn. The 
conclusion is that the presence of the fermion determinant completely alters 
the vacuum structure of the theory. The question we wish to address is how 
we can understand this based on the spectrum of the QCD Dirac operator. For 
simplicity we only consider the case of Nf = 1 and = 0. 

The average spectral density in full QCD is defined by 

/""(x, y, fi) = S\x + iy- Xk) det{D + m + njo))- (9.1) 
k 

The low-energy limit of the generating function for the spectral density can 
again be written as a r-function [26]. The spectral density is then obtained 
from the replica Umit of the corresponding Toda lattice equation. The result is 

full. ^ y3^x2 + 2/2)s4 F(a;2+y2)22 
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It was first obtained from the random matrix model [27] using the method 
of complex orthogonal polynomials developed in [89, 90]. To appreciate this 
result, we considers its asymptotic expansion for F ^ oo. For m = and 
x > we can derive the asymptotic result for the difference 

p^"-(ar, y, p) - p'-\x, y, p) ~ W^'^^e^*^^. (9.3) 
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The behavior near the extremum at (.x, y) = (4^^F^/3S, 0) is given by 

y, m) - P^''"(a;, y, fi) ~ ei'^'-^'^ei^^^^. (9.4) 

We have oscillations on the scale 1/TjV with an amplitude that diverges expo- 
nentially with the volume [26]. In the thermodynamic limit, these oscillations 
are visible in a domain where the real part of the exponent in (9.3) is positive. 
This region is given by intersection of the inside of the ellipses 

3{x±^^^'F^Ef + y^=^-^^^^F^J:^ (9.5) 

and the strip \x\ < 2F^/l(^/E. At the mean field level this can be reinterpreted 
as a region where Kaon condensation takes place [94, 95]. 

10. Conclusions 

The existence of two formulations of QCD at low energy, first as a mi- 
croscopic theory of quarks and gluons and second as an effective theory of 
weakly interacting Goldstone bosons, imposes powerful constraints on either 
of the theories. The effective theory is completely determined by the symme- 
tries of the microscopic theory, and the mass dependence of the effective theory 
imposes sum rules on the inverse Dirac eigenvalues. In particular this means 
that any theory with the same symmetry breaking pattern and a mass gap will 
be subject to the same constraints. The simplest microscopic theory is chiral 
Random Matrix Theory. 

However, more can be done than constraining the inverse Dirac eigenvalues 
by sum rules. The key observation is that the generating function for the resol- 
vent amounts to a partition function with additional flavors with a mass z equal 
to the value for which the resolvent is calculated. Again we have a microscopic 
theory and an effective theory with the same low energy limit. Because z is 
a free parameter, it can always be chosen such that the Compton wavelength 
of the corresponding Goldstone boson is much larger than the size of the box. 
In this region the z dependence of the partition function is determined by the 
mass term of the chiral Lagrangian which is a simple matrix integral. 

To obtain the Dirac spectrum we have to quench the determinant corre- 
sponding to z. This can be done in two ways: by the replica trick or by the 
supersymmetric method. Although, the supersymmetric method is straightfor- 
ward, the naive replica trick is technically somewhat simpler. The problem is 
that the rmive replica trick does not give correct nonperturbative results. One 
way out is if the dependence on the number of replicas n is known as an an- 
alytical function of n around n = 0. It was shown by Kanzieper that the 
n-dependence can be obtained from the solution of the Painleve equation. The 
other way out, which has been advocated in these lectures, is if the partition 
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function for n = is related by a recursion relation to partition functions 
with a nonzero integer number of flavors. The replica limit of this Toda lattice 
equation gives us nonperturbative correlation functions. This is an efficient 
formulation of the problem. The sttucture of the final answer already has the 
factorized structure of the Toda lattice formulation. We could also say that the 
supersymmetric partition function connects two semi-infinite hierarchies. 

New results with the Toda lattice method were obtained for QCD at nonzero 
chemical potential. In this case the low-energy effective partition functions are 
also related by the Toda lattice equation. This made it possible to express the 
microscopic spectral density as the product of the partition function with one 
fermionic flavor and the partition function with one bosonic flavor. This result 
has later been reproduced by RMT with the method of orthogonal polynomials. 

More surprisingly the Toda lattice method also gives the correct result for 
QCD at nonzero chemical potential with dynamical fermions. Because of the 
phase of the fermion determinant, a breakdown of this method for this case 
would not have been a surprise. However, the concept of integrability that also 
reigns this case, is so powerful that replica limit can be taken in exactly the 
same way as in the quenched case. The result for the spectral density shows 
oscillations on the scale of 1/F and an amplitude that diverges exponentially 
with V. This structure is necessary to obtain a nonzero chiral condensate in 
the chiral limit. 

The Toda lattice method has been applied to quite a few cases in the sym- 
metry class /3 = 2 (see [24, 23, 92, 25, 26, 93]). Our conjecture is tiiat all 
microscopic correlation functions in this class can be obtained from the replica 
hmit of a Toda lattice equation. A much tougher problem is analysis of the 
replica limit for the other Dyson classes. We are not aware of any progress on 
this problem and encourage the reader to confront this challenge. 
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